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We consider the transition of the Nagel-Schreckenberg traffic flow model from the free flow regime 
to the jammed regime. We examine the inhomogeneous character of the system by introducing a 
new method of analysis which is based on the local density distribution. We investigated the 
characteristic fluctuations in the steady state and present the phase diagram of the system. 
PACS numbers: 89.40.-Fk,05.40.-Fj,05.60.-Hw 
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I. INTRODUCTION 

Over the past few years much attention has been de- 
voted to the study of trafhc flow. Since the seminal 
work of Lighthill and Whitham in the middle of the 50 's 
[01 many attempts have been made to construct more 
and more sophisticated models which incorporate various 
henomena occurring in real traffic (for an overview see 
). Recently, a new class of models, based on the idea 
of cellular automata, has been proven to describe traffic 
dynamics in a very efficient way [^ . Especially the tran- 
sition from free flow to jammed traffic with increasing car 
density could be investigated very accurately. Neverthe- 
less, besides various indications [^ , no unique description 
for a dynamical transition could be found. Furthermore, 
no satisfying order parameter could be defined so far. 
In this article we introduce a new method of analysis 
which allows us to identify the different phases of the 
system and to describe the phase transition in detail, 
i.e., defining an order parameter, considering the fluctu- 
ations which drive the transition, and determining the 
phase diagram. 

We consider a one-dimensional cellular automaton of 
linear size L and N particles. Each particle is as- 
sociated the integer values Vi G {0,1,2, ...,Vm,ax} and 
di £ {0, 1, 2, 3, ...}, representing the velocity and the dis- 
tance to the next forward particle [^. For each particle, 
the following four update steps representing the acceler- 
ation, the slowing down, the noise, and the motion of 
the particles are done in parallel: (1) if < di then 
Vi Miii{vi + l,Vr,^ax}, (2) if Vi > di then Vi di, 
(3) with probability P Vi ^ Ma,x{vi — 1,0}, and (4) 
ri Ti + Vi, where denotes the position of the i-th 
particle. 



II. SIMULATIONS AND RESULTS 

Figure ^ shows a space-time plot of the system. Each 
dot corresponds to a particle at a given time step. The 
global density pg = N/L exceeds the critical density and 
jams occur. Traffic jams are characterized by a high lo- 
cal density of the particles and by a backward movement 
of shock waves [|l|. One can see from Fig.|] that in the 
jammed regime the system is inhomogeneous, i.e., traf- 
fic jams with a high local density and free flow regions 
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FIG. 1. Space-time plot for Vmax = 5, P = i, and pg > pc. 
Note the separation of the system in high and low density 
regions. 



with a low local density coexist. In order to investigate 
this transition one has to take this inhomogeneity into 
account. 

Traditionally one determines the so-called fundamen- 
tal diagram, i.e., the diagram of the flow vs the density. 
The global flow is given by, ^ — pg (v) , where (v) 
denotes the averaged velocity of the particles. Due to 
the stochastic behavior of the dynamics for < P < 1 
the system behaves independently of the initial condi- 
tions after a certain transient regime. In this limit one 
can interpret (...) as a time or ensemble average. For 
P = and P = 1 the dynamics is deterministic and the 
behavior depends strongly on the initial conditions. In 
any case, this non-local measurements are not sensitive 
to the inhomogeneous character of the system, i.e., the 
information about the two different coexisting phases is 
lost. In the following we apply a new method of analysis 
which is based on the measurement of the local density 
distribution p{p). The local density p is measured on a 
section of the system of size 5 according to 
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FIG. 2. The local density distributioijjp(p) for various val- 
ues of the global density, Vmax — 5, P — ^ and S = 256. The 
dashed line corresponds to the characteristic density of the 
free flow phase. 



FIG. 3. The maximum of the local density distribution as 
a function of the global density pg. The inset shows that 
the transition point pc does not depend on the value of the 
parameter S. 



Of course we have checked that the main results are not 
affected by the value of S, provided that delta is signifi- 
cantly smaller than the system size L in order to measure 
the local properties. In order to reflect the behavior of 
the low density regime 6 should be significantly larger 
than a certain length scale Aq which corresponds to the 
characteristic length scale of the density fluctuations in 
the free flow phase (see below). For any parameter set 
{vmax, P} the local density p fluctuates around the value 
of the global density pg and the probability distribution 
of the local density p{p) contains all informations needed 
to describe the transition. 

The local density distribution p{p) is plotted for vari- 
ous values of the global density Pg in Fig. |^. In the case 
of small values of pg, see Fig. ^ja, the particles can be 
considered as independent (see below) and the local den- 
sity distribution is simply Gaussian with the mean val- 
ues Pg and a width which scales with ^/6. Increasing the 
global density, jams occur and the distribution displays 
two different peaks (Fig. |^). The first peak corresponds 
to the density of free particles and in the phase coexis- 
tence regime the position of this peak does not depend 
on the global density (see the dashed lines in Fig. ||). 
The second peak is located at larger densities and char- 
acterizes the jammed phase. With increasing density the 
second peak occurs in the vicinity of the critical density 
Pc (Fig. ||b) and grows further (Fig. ||c) until it dom- 
inates the distribution in the sense that the first peak 
disappears (Fig. ||d). The two peak structure of the local 
density distribution clearly reflects the coexistence of the 
free flow and jammed phase above the critical value Pc- 
In the following we show that the behavior of the peaks 
leads to a determination of pc- 

One expects that in a homogeneous system the local 
density distribution displays one peak and is symmetric 
around the global density, i.e., p{pmax) = Pg- This can- 



not be the case in an inhomogeneous system where the 
local density distribution displays two peaks correspond- 
ing to two different coexisting phases. In Fig. ^ we plot 
the position of the maximum of the local density distri- 
bution pipmax) as a function of the global density pg. 
One clearly sees the transition point pc where the posi- 
tion of the maximum becomes independent of the global 
density. The inset of Fig. ^ shows that the determina- 
tion of the transition point does not depend on the special 
value of S. Only the point where the second peak exceeds 
the first peak depends on the measurement parameter 6. 
With increasing 6 this point tends to smaller values of 
Pg because with increasing S the measurement starts to 
average over the two different phases. From these mea- 
surements we conclude that the phase transition of the 
Nagel-Schreckenberg model is a transition from a homo- 
geneous regime (free flow phase) to an inhomogeneous 
regime which is characterized by a coexistence of two 
phases (free flow traffic and jammed traffic). 

In order to describe the spatial decomposition of the 
coexisting phases we measured the steady state structure 
factor (I 



Sik) 



(2) 



where 77 (r) = 1 if the lattice site r is occupied and 
ri{r) — otherwise. In Fig. ^ we plot the structure fac- 
tor S{k) for the same values of the global density as in 
Fig. i.e., below, in the vicinity, above and far away of 
the transition point. It is remarkable that S{k) exhibits a 
maximum for all considered values of the global density 
at ko « 0.72 (dashed lines in Fig. |^). This value cor- 
respondence to the characteristic wave length Aq = |^ 
of the density fluctuations in the free flow phase. The 
steady state structure factor is related to the Fourier 
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FIG. 4. The structure factor S(fc) for = |, Vmai. = 5 and 
for various values of the density p. The dashed line marks the 
characteristic wavelength Ao of the free flow phase. 



transform of the real space density-density correlation 
function. The "wave length Aq corresponds to a maximum 
of the correlation function, i.e., Aq describes the most 
likely distance of t"wo particles in the free flow phase. For 
lo"W" densities the structure factor is almost independent 
of the density and displays a minimum for small k values 
indicating the lack of long-range correlations. Crossing 
the transition point the smallest mode S{k — ^) in- 
creases quickly. This suggests that the jammed phase is 
characterized by long-range correlations which decay in 
the limit Pg ^ pc algebraically as one can see from the 
log-log plot in Fig. ||d. 

We already mentioned that fco characterizes the den- 
sity fluctuations in the free flow phase. In Fig. || we 
plot fco as a function of Vmax- Except for the case 
Vraax — 1 the maximum velocity and fco obey the relation 



fco (w„ 



1) 



const. The fact that ti„ 



1 violates 



this equation does not surprise. It is already known that 
the physics for Vmax > 2 is distinctly different from the 
case Vmax = 1 @- For instance, for Vmax = 1 the fun- 
damental diagram is symmetric around its maximum at 
Pg — 0.5 independent of the noise parameter P. Whereas, 
the position of the maximum depends on P for Vmax > 2 
and no symmetry occurs. Another example is that for 
Vmax > 2 jams are allowed to branch (see Fig. |i), unlike 
jams for Vmax = 10- The qualitative different behavior 
for Vmax = 1 is caused by a particle-hole symmetry which 
is lost for larger values of the maximum velocity . 

Changing from momentum space to real space the 
characteristic wavelength Ao of the density fluctuations 
in the free flow phase are given by 



Ao 



27r 
fc'o 



27r. 



const 



(3) 



On the other hand the average distance d of the parti- 
cles is given by the inverse density p~^. In the free flow 
phase the average distance is larger than the wavelength 



FIG. 5. The characteristic inverse wavelength fco for the 
maximum velocities Vmax G {1, 2, 3, 4, 5, 8} for P = |. 



Ao, Ao ^ d — 1/ Pg, i.e., the cars can be considered 
as independent particles. With increasing density this 
behavior changes when the average distance d is compa- 
rable to the wavelength Ao ~ d. The critical density is 
related to the characteristic wavelength 



Pc 



1 

A^ 



1 



1 



(4) 



1 is 



The fact that the critical density scales with Vm, 
already known for the deterministic case P = 

Up to now we only considered the case P — \- The 
phase diagram in Fig. || shows the P dependence of the 
transition density pc- /denotes the free flow phase and 
f+i corresponds to the coexistence region where the sys- 
tem separates in the free flow and jammed phase. The 
dashed line displays the P dependence of the maximum 
flow obtained from an analysis of the fundamental di- 
agram 1^. The critical densities pc, where the phase 
transition takes place, are lower than the density values 
of the maximum flow. Measurements of the relaxation 
time, which is expected to diverge at a transition point 
Q, confirm this result (see Fig. |^). But one has 
to mention that the determination of the critical density 
via relaxation times leads in the coexistence regime f+j to 
unphysical results, in the sense that the relaxation time 
becomes negative [ pl)| . 

Motivated by real traffic fiow we fixed in our measure- 
ments the noise parameter P and varied the global den- 
sity Pg, i.e., we crossed the critical line in the phase dia- 
gram parallel to the vertical axis (see Fig. ||). In theoret- 
ical investigations however, it is more convenient and re- 
vealing to consider the crossing of the critical line parallel 
to the horizontal axis, i.e., increasing the noise parame- 
ter P by a fixed density. With growing P density fiuc- 
tuations with a characteristic wavelength Ao occur. This 
wavelength Ao grows with increasing P but no phase sep- 
aration takes place until it exceeds at the transition line 
a critical wavelength Ac. Then the system separates into 
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FIG. 6. The phase diagram of the Nagel-Schreckenberg 
model. Note that in the non-deterministic region < P < 1 
the density of the maximum flow exceeds the density of the 
transition point. 
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FIG. 7. The number of particles in the high density phase 
Nj normalized by the total number of particles A''. Nj/N 
could serve as an order parameter and the critical exponent 
is /? = 1 (see solid line). 



the two coexisting phases and the amount of particles 
which belong to the jammed phase, Nj, could serve as 
an order parameter. As mentioned above the local den- 
sity distribution displays two peaks in the heterogeneous 
phase and the area under each peak is proportional to Nf 
and Nj, respectively. Figure ^ displays Nj normalized 
by the total number of particles {N — Nf + Nj). Ap- 
proaching the transition point p^, Nj arises linear, i.e., it 
obeys the equation Nj ~ (P — Pc)'^ , where the critical 
exponent is given by P — 1. The continuously behav- 
ior of Nj and the diverging relaxation time ||l^ suggest 
that the transition of the traffic model could describe as 
a phase transition of second order. A detailed analysis 
of the order parameter and the order parameter fluctu- 
ations, including a finite-size analysis, requests further 
investigations. 

III. CONCLUSIONS 
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In conclusion we have studied numerically the Nagel- 
Schreckenberg traffic flow model using a local density 
analysis. Crossing the critical line of the system a phase 
transition takes place from a homogeneous regime (free 
flow phase) to an inhomogeneous regime which is char- 
acterized by a coexistence of two phases (free flow traffic 
and jammed traffic). The decomposition in the phase 
coexistence regime is driven by density fluctuations, pro- 
vided they exceed a critical wavelength Ac. The amount 
of particles in the jammed phase could serve as an or- 
der parameter which arises linear at the transition point, 
suggesting that the transition is of second order. 



